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Abstract 

The chloroplasts of land plants contain internal membrane systems, the thylakoids, which are arranged in stacks 
called grana. Because grana have not been found in Cyanobacteria, the evolutionary origin of genes controlling the 
structural and functional diversification of thylakoidal membranes in land plants remains unclear. The angulata10-1 
(anu10-1) mutant, which exhibits pale-green rosettes, reduced growth, and deficient leaf lateral expansion, resulting 
in the presence of prominent marginal teeth, was isolated. Palisade cells in anu10-1 are larger and less packed than 
in the wild type, giving rise to large intercellular spaces. The ANU10 gene encodes a protein of unknown function 
that localizes to both chloroplasts and amyloplasts. In chloroplasts, ANU10 associates with thylakoidal membranes. 
Mutant anu10-1 chloroplasts accumulate H 2 0 2 , and have reduced levels of chlorophyll and carotenoids. Moreover, 
these chloroplasts are small and abnormally shaped, thylakoidal membranes are less abundant, and their grana are 
absent due to impaired thylakoid stacking in the anu10-1 mutant. Because the trimeric light-harvesting complex II 
(LHCII) has been reported to be required for thylakoid stacking, its levels were determined in anu10-1 thylakoids and 
they were found to be reduced. Together, the data point to a requirement for ANU10 for chloroplast and mesophyll 
development. 

Keywords: Arabidopsis thaliana, chloroplast, grana, LHCII trimers, mesophyll development, thylakoid biogenesis, thylakoid 
stacking. 



Introduction 

In land plants, mutants with defective pigment content have 
allowed the identification of numerous nuclear genes that 
are crucial for organelle division and other plastid-specific 
processes (Lopez- Juez and Pyke, 2005; Hricova et ah, 2006). 
Many of these genes encode plant-specific proteins whose 
closest homologues are found in Cyanobacteria, as expected 
if they have been acquired by the nuclear genome in a hori- 
zontal transfer event, originating in the genome of the cyano- 
bacterial ancestor of plastids. As an example, the nuclear 
genomes of land plants encode homologues of Filamenting 
temperature-sensitive Z (FtsZ) (Osteryoung and Vierling, 
1995; Osteryoung et ah, 1998), Minicell D (MinD) (Colletti 
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et ah, 2000), MinE (Itoh et al, 2001), and many other com- 
ponents of the prokaryotic cell division apparatus. 

The chloroplasts of land plants contain internal mem- 
brane systems, the thylakoids, which are arranged in stacks 
called grana. Grana thylakoids form cylindrical stacks that 
are connected to adjacent grana by non-stacked stroma thy- 
lakoids, which form right-handed helices around the grana 
(Austin and Staehelin, 2011). The thylakoidal membranes 
harbour the photosynthetic protein complexes, photosys- 
tems I and II (PSI and PSII), and other protein complexes, 
including ATP synthase, cytochrome bff, and the light- 
harvesting complexes I and II (LHCI and LHCII). Grana 
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underlie a differential specialization of thylakoid membranes 
in the lateral dimension: the PSII and LHCII complexes are 
more abundant in grana, and the PSI and ATP synthase 
are preferentially found in the adjacent stroma thylakoids 
(Albertsson, 2001; Dekker and Boekema, 2005). Because 
grana have not been found in Cyanobacteria, the evolution- 
ary origin of genes controlling the structural and functional 
diversification of thylakoidal membranes in land plants 
remains unclear (Mullineaux, 2005). Mutants with abnormal 
thylakoid stacking and loss of typical grana have long been 
known in species such as maize (Bachmann et al, 1969), 
barley (Nielsen et al, 1979), tobacco (Archer and Bonnett, 
1987), and Arabidopsis thaliana (hereafter, Arabidopsis) 
(Reiter et al. , 1994). However, the molecular basis of the phe- 
notypes observed in these species has only been determined 
in a few cases. Mutations in the Arabidopsis CHLORATA-42 
(CH-42) gene (Koncz et al, 1990), which encodes one of 
the three subunits of the chloroplast magnesium chelatase 
complex, significantly reduce thylakoid stacking, suggesting 
a link between chlorophyll biosynthesis and grana forma- 
tion (Apchelimov et al, 2007). Another Arabidopsis gene, 
GRANA-DEFICIENT CHLOROPLAST1 {GDC1), encodes 
an ankyrin-domain protein that is essential for grana for- 
mation. GDC1 is required for the assembly of the trimeric 
forms of the LHCII complex, which are barely detected in 
the gdcl-3 mutant (Cui et al. , 201 1). Current models empha- 
size the role of LHCII trimers in the formation of grana and 
the stacking of thylakoidal membranes (Mullet and Arntzen, 
1980; Day et al, 1984; Garab and Mustardy, 1999; Allen 
and Forsberg, 2001), mainly through electrostatic interac- 
tions between the positively charged N-terminal domains of 
LHCII trimers (stroma-exposed) and the negatively charged 
surface of LHCII trimers in adjacent thylakoidal membranes 
(Steinback et al, 1979; Standfuss et al, 2005). In addition, 
members of the CURVATURE THYLAKOID 1 (CURT1) 
family of proteins have recently been shown to be enriched 
at the margin of grana, where they are thought to promote 
the curvature of thylakoidal membranes in Arabidopsis. 
Interestingly, the Arabidopsis CURT1A protein can replace 
the function of a distant cyanobacterial orthologue, show- 
ing that at least some proteins with crucial roles in thyla- 
koid architecture are evolutionarily conserved, even though 
Cyanobacteria lack grana (Armbruster et al, 2013). 

Here it is shown that the ANGULATA10 (ANU10) gene 
of Arabidopsis encodes a novel plastid-localized protein 
that is conserved among land plants, and the phenotypic 
and molecular characterization of loss-of-function anulO 
mutants is reported. It is demonstrated that the ANU10 
protein localizes to plastids, including amyloplasts and chlo- 
roplasts, where it is associated with thylakoidal membranes. 
Thylakoid biogenesis is seriously impaired in anul0-l chlo- 
roplasts. Supporting the implication of LHCII trimers in 
grana formation, the anul0-l mutant contains reduced levels 
of LHCII trimers and shows defective thylakoid stacking, 
therefore lacking typical grana. A relationship between plas- 
tid integrity and leaf development can be inferred from the 
larger, sparsely packed cells observed in the palisade meso- 
phyll of the anul0-l mutant. 



Materials and methods 

Plant material and growth conditions 

Arabidopsis thaliana (L.) Heynh. wild-type accessions Columbia-0 
(Col-0) and Landsberg erecta (Ler), as well as the T-DNA insertion 
lines SAIL_708_F05 (N831342) and SAIL_659_F07 (N828696), 
were obtained from the Nottingham Arabidopsis Stock Centre 
(NASC). The anul0-l mutant was isolated after ethylmethane sul- 
phonate (EMS (-induced mutagenesis as previously described (Berna 
et al, 1999). All plants in this work were grown on Murashige and 
Skoog agar medium (2.15g l" 1 ), at 20 ±1 °C and 60-70% relative 
humidity under continuous fluorescent light (~90 umol m 2 s~') as 
previously described (Ponce et al, 1998). For the reactive oxygen 
species (ROS) production assay, low intensity (~55 umol mT 2 s~') 
and moderately high intensity (~180 umol nT 2 s~') light conditions 
were additionally used. Crosses and allelism tests were performed as 
reported in Berna et al. (1999). 



Positional cloning and molecular characterization of anul 0 
alleles 

Low-resolution mapping of the anul0-l mutation was performed 
as described in Ponce et al. (2006). For fine mapping, the nga392, 
S0392, and cer479911 insertion/deletion polymorphisms from 
Monsanto (http://www.arabidopsis.org/browse/Cereon) and the 
F3M18 and F1K23 single nucleotide polymorphisms (SNPs) from 
the 1001 genomes project database (Weigel and Mott, 2009) were 
used. To find the anul 0-1 mutation, a 39 18 bp fragment encom- 
passing the entire transcription unit of Atlg28530 was PCR ampli- 
fied from Ler and anul0-l genomic DNA, and sequenced on an 
ABI PRISM 3130x1 Genetic Analyser (Applied Biosystems). All 
the primers used for the cloning and sequencing of Atlg28530 
are listed in Supplementary Table SI available at JXB online. 
The T-DNA insertions in N831342 and N828696 lines were 
confirmed by PCR using primers recommended by the 'T-DNA 
Primer Design' tool (http://signal.salk.edu/tdnaprimers.2.html; 
Supplementary Table SI). 



Bioinformatic analyses 

Four cDNA (BT005784.1, BT008607.1, BX816836.1, and 
AK228671.1) and nine expressed sequence tag (EST; AU228187.1, 
AU237163.1, BP599159.1, BP606540.1, BP662266.1, BP807381.1, 
BP846942. 1 , ES01 7536. 1 , and EL9842 16.1) sequences from GenBank 
were used to assemble the transcriptional unit of Atlg28530 using 
CAP3 (http://pbil.univ-lyonl.fr/cap3.php) (Huang and Madan, 
1999). Subcellular localization was predicted with TargetP (http:// 
www.cbs.dtu.dk/services/TargetP/) (Emanuelsson et al, 2000) and 
Multiloc2 (http://abi.inf.uni-tuebingen.de/Services/MultiLoc2) 
(Blum et al, 2009). Chloroplast transit peptide sequences and 
transmembrane domains were predicted with ChloroP 1.1 (http:// 
www.cbs.dtu.dk/services/ChloroP/) (Emanuelsson et al, 1999) and 
SOSUI (http://bp.nuap.nagoya-u.ac.jp/sosui/) (Hirokawa et al, 
1998), respectively. 

To identify ANU10 homologues, BLASTP searches (Altschul 
et al. , 1997) were carried out at the NCBI server using a word size of 
2 and default values for all other parameters. Full-length sequences 
were selected based on a BLAST E- value cut-off of 3 x 10 -6 , and 
were subsequently aligned using the consistency-based method 
implemented in T-Coffee (Notredame et al, 2000). Alignments 
were refined with MUSCLE (http://www.ebi.ac.uk/Tools/msa/mus- 
cle/) (Edgar, 2004) and shaded with BOXSHADE3.21 (http://www. 
ch.embnet.org/software/BOX_form.html). Identity percentages 
were calculated using the Sequence Manipulation Suite (http://www. 
bioinformatics.org/sms2/index.html) (Stothard, 2000). Phylogenetic 
trees were obtained using MEGA5 (Tamura et al, 2011). Searches 
for distant homologues of ANU10 were carried out with HMMER 
(http://hmmer.janelia.org/) (Finn et al. , 201 1). 
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RNA isolation, cDNA synthesis, and qRT-PCR 

Total RNA from Col-0 rosettes [collected 21 d after stratification (das)] 
was extracted using TRI Reagent (Sigma), and DNA was removed 
using the TURBO DNA-free Kit (Invitrogen). First-strand cDNA 
was synthesized using random hexamers and the Maxima Reverse 
Transcriptase system (Fermentas). For quantification of the expres- 
sion of nuclear and plastid genes, the primers listed on Supplementary 
Table SI available at JXB online were used. The 18S rRNA gene was 
used as an internal control (Yamauchi et al. , 2004). Three different 
biological replicates and triplicate reactions were used. Amplification 
reactions were prepared in a volume of 20 ul by adding 7.5 ul of 
Maxima SYBR Green/ROX qPCR Master Mix (Fermentas), 3 ul of 
the corresponding primer pair (2.5 uM each), and 1 ul of cDNA tem- 
plate. Relative quantification of gene expression data was performed 
using the comparative C T method (Schmittgen and Livak, 2008) on a 
Step One Plus System (Applied Biosystems). 

Gene constructs and plant transformation 

To make the 35S pm :ANUJ0 and 35S pm :ANU10:GFP (green fluo- 
rescent protein) transgenes, the full-length coding sequence of 
ANUI0 was amplified from Col-0 cDNA using Phusion polymerase 
(Thermo Scientific) with the ANUlOcds-F and ANUlOcds-R prim- 
ers (Supplementary Table SI at JXB online). The amplification prod- 
uct was cloned into the pENTR/D-TOPO entry vector (Invitrogen) 
and transferred into the pMDC32 and pMDC83 destination vec- 
tors, which include a dual 35S promoter (Curtis and Grossniklaus, 
2003). For the ANU10 pm :GUS ^-glucuronidase) transgene, a 1.5-kb 
fragment encompassing the intergenic region between Atlg28530 
and Atlg28540, including the Atlg28530 5'-untranslated region 
(UTR) and Atlg28540 3'-UTR, was amplified using Col-0 genomic 
DNA as the template. The resulting product was cloned into the 
pGEM-T Easy221 vector (kindly provided by B. Scheres), and 
then transferred into the pMDC163 destination vector (Curtis and 
Grossniklaus, 2003). Alternatively, the ANU10 promoter region was 
fused with the ANU10 full-length cDNA in a two-template PCR 
of overlapping products, using the primers ANUlOcds-pro and 
ANUlOpro-cds (Supplementary Table SI). The resulting fusion 
product was cloned into the pGEM-T Easy221 vector and trans- 
ferred into pMDCl 1 1 (Curtis and Grossniklaus, 2003) to obtain the 
ANU10 pro :ANU10:GFP construct. All the constructs were trans- 
formed into Agrobacterium tumefaciens LBA4404. Col-0, Ler, and 
anul0-l plants were transformed by the floral dip method (Clough 
and Bent, 1998). Tj transgenic plants were selected on plates sup- 
plemented with 15 ug mr 1 hygromycin B (Invitrogen). 

Microscopy, histology, and morphometry 

For transmission electron microscopy, mutant, wild-type, and 
transgenic plants were harvested 16 das. Leaf tissue excluding the 
primary vein and leaf margin was excised and fixed in McDowell's 
solution (McDowell and Trump, 1976), and prepared as previously 
described (Hricova et al., 2006). Samples were visualized at 80 kV 
using a JEOL 1011 transmission electron microscope equipped with 
a Gatan 792 BioScan digital camera. Grana diameter and height 
measurements were obtained from transmission electron micro- 
graphs using the ImageJ software. Leaf tissues were imaged using a 
Leica TCS SPE confocal microscope. Root tissues were imaged using 
a Nikon CI confocal microscope. Transverse sections of leaves were 
obtained as described by Serrano-Cartagena et al. (2000), embed- 
ding the tissue in Technovit7100 resin and obtaining 10 urn sections. 
Rosette area measurement and morphometric analysis of palisade 
cells and transverse sections of leaves were performed as described 
previously (Perez-Perez et al, 2011; Ferrandez-Ayela et al, 2013). 

Dry weight and pigment determination 

For dry weight measurement, eight plants of each genotype were 
oven-dried for 48 h at 55 °C. For determination of chlorophylls and 



carotenoids, four independent samples of lOOmg of fresh leaves 
from rosettes collected 16 das were pooled, frozen in liquid N 2 , and 
homogenized with 3.5ml of cold 80% acetone. The samples were 
centrifuged for 5min at 5000 rpm and the pigment concentration in 
the supernatant was spectrophotometrically determined as previ- 
ously described (Wellburn, 1994). 



ROS, GUS, and Lugol staining 

To determine ROS accumulation, a minimum of six first-node 
leaves from each genotype were excised and incubated overnight 
in 1 mg mL 1 3,3'-diaminobenzidine (DAB; Sigma-Aldrich), under 
vacuum and in the dark. Leaves were cleared by boiling in acetic 
acid:glycerol:ethanol (1:1:3 v/v/v) for 5min, and then incubated in 
96% (v/v) ethanol until chlorophyll was bleached. GUS staining 
of plant tissues and Lugol staining of root tips were performed as 
described in Robles et al. (2010) and Willemsen et al (1998), respec- 
tively. Samples were visualized with a Nikon CI microscope. 



Western blot analysis of chloroplast proteins and blue native 
PAGE of thylakoid membranes 

Chloroplasts from rosette leaves collected 16 das were isolated as 
described by Grabsztunowicz and Jackowski (2012), omitting the 
Percoll step gradient. Stroma and thylakoid membranes from iso- 
lated chloroplasts were separated as previously reported ( Armbruster 
et al. , 2010). A volume of isolated chloroplasts and thylakoids equiv- 
alent to 20 ug of chlorophyll, as well as 20 ug of proteins from the 
stroma fraction, including a lane with an EZ-Run Pre-stained Rec 
Protein Ladder (Fisher BioReagents), were resolved by 10% SDS- 
PAGE, blotted on nitrocellulose membranes (Amersham Hybond 
ECL, RPN203D; GE Healthcare), and subjected to immunoblot 
analysis with specific antibodies (G1544, Sigma; AS03 037-10 and 
AS05 092, Agrisera AB). Isolation of thylakoids from rosettes col- 
lected 16 das and blue native PAGE were carried out as reported in 
Perez-Perez et al. (2013). 



Results 

Leaf and whole-plant defects in anu10-1 plants 

In a large-scale screen for EMS-induced mutants with abnor- 
mal leaf shape, 18 mutants with pale-green leaves and dentate 
margins were previously identified. The causative mutations 
were dubbed angulatci (anu), and complementation tests 
showed that they damage 12 different genes (ANU1-ANU12) 
(Berna et al, 1999). 

The only mutant allele of the ANU10 gene, anul0-l, 
causes reduced rosette and stem growth (Fig. 1A, B, J). The 
dry weight and projected area of anul0-l rosettes were sig- 
nificantly lower than in wild-type Ler plants throughout 
the study (Fig. IK; Supplementary Fig. SI at JXB online). 
Rosettes of the anul0-l mutant harvested 21 das showed a dry 
weight of 3.02±0.16mg and an area of 147.72 ±32.68 mm 2 , 
while the corresponding values for Ler were 6.02 ± 0.91 mg 
and 701.48 ± 56.88mm 2 , respectively. Adult anul0-l plants 
also exhibited a significant reduction in shoot length: 
28.47±2.69cm in her, but only 18.22±3.55cm in anulO- 
1 42 das (Fig. 1J; Supplementary Fig. SI). The number of 
secondary stems was lower in anul0-l (4.67 ±0.52) than in 
Ler (9.25 ±0.71) determined 35 das, although this difference 
disappeared over time (Supplementary Fig. SI). The leaves 
of anul0-l were pale green and their margins had prominent 
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Fig. 1. Phenotypic characterization and rescue of the anu10-1 mutant. (A-C) Rosettes, (D-F) first-node leaves, and (G-l) bright-field micrographs of 
the subepidermal layer of palisade mesophyll cells from (A, D, G) the Ler wild type, (B, E, H) the anu10-1 mutant, and (C, F, I) a transgenic anuW-1 
35Spro:ANU10 plant. (J) Adult plants. Pictures were taken (A-l) 16 and (J) 42 das. Scale bars indicate (A-F) 2 mm, (G-l) 30 \im, and (J) 5cm. (K) Dry 
weight and (L) pigment content in Ler, anuW-1 , and transgenic anul0-1 35Spro:ANU10 plants. Error bars indicate standard deviations. Asterisks 
indicate values significantly different from Ler in a Mann-Whitney U-test [(K) P<0.01 , n=8 and (L) P<0.05, n=4]. 



teeth (Fig. IE). As seen in paradermal sections, the palisade 
mesophyll cells were irregular in size, with larger cells in 
anulO-1 than in the wild type (Fig. 1G, H). By measuring the 
size of individual cells, a shift in the distribution of palisade 
cell size towards greater values (Fig. 2A), as well as the pres- 
ence of large intercellular spaces in the palisade mesophyll 
of anul0-l leaves (Fig. 2B, C, E, F) was detected. Similar 
defects were observed in transverse sections of first-node 
leaves (Fig. 2H, I). In these sections, mesophyll cells filled a 
significantly smaller percentage of the section area, which 
was matched by a significant increase in the area occupied by 
air spaces in anul0-l leaves (Fig. 2K). 

To gain insight into the physiological basis of the pale- 
green phenotype of anul0-l, pigment levels were measured 
in the mutant. In line with the observed pale-green pheno- 
type (Fig. IB, E), a reduction in the levels of chlorophyll a, 



chlorophyll b, and carotenoids was detected (Fig. 1 L). Because 
carotenoids play an important role in photoprotection, the 
presence of ROS, visualized by staining with DAB, was also 
tested. Darker staining was observed in amil0-l leaves rela- 
tive to Ler (Supplementary Fig. S2 available at JXB online), 
which indicates the presence of increased levels of H 2 0 2 . In 
the mutant, the accumulation of H 2 0 2 was positively corre- 
lated with light intensity (Supplementary Fig. S2). 

Positional cloning of ANU 10 

The anul0-l mutation was previously mapped to chromo- 
some 1 (Robles and Micol, 2001). To understand further the 
molecular basis of the phenotype of anul0-l, a positional 
cloning approach was undertaken to identify the causal gene. 
First a 72 kb candidate interval (Fig. 3A) flanked by two SNP 
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markers, F3M18 and F1K23 (Supplementary Table SI avail- 
able at JXB online), was denned. This interval encompassed 
24 candidate genes. Because the phenotype of anul0-l sug- 
gests a defect in a chloroplast-related function (Fig. IB, E), 
the Atlg28530 gene, which is the only gene in this interval 
predicted to encode a chloroplast-localized protein, was the 
focus of further study. The Atlg28530 gene was sequenced 
in mutant and wild-type plants, and a G— >A transition was 
found in its coding region only in anul0-l mutant plants 
(Fig. 3B). Expression of the Atlg28530 gene is supported by 
four cDNA and nine EST sequences deposited in GenBank 
(see the Materials and methods). In order to determine the 
intron-exon structure of the Atlg28530 gene, these sequences 
were assembled using CAP3 and the consensus mRNA 
sequence was aligned to the sequence of Arabidopsis chro- 
mosome 1 (NC_003070). The consensus mRNA sequence 
has eight exons and its longest open reading frame encodes 
a protein with 614 amino acids and a molecular mass of 
68.46 kDa. The anul0-l mutation introduces a stop codon 
(Trp317— >Stop) in the third exon of the gene (Fig. 3B), and 
is predicted to truncate the protein prematurely (from 614 
amino acids in her to only 316 in anul0-l). 

To identify additional alleles of the ANU10 gene, two 
publicly available lines, SAIL_708_F05 and SAIL_659_F07, 
which carry T-DNA insertions in the coding sequence of 
Atlg28530 were characterized. When homozygous for the 
T-DNA insertion, these two lines displayed phenotypes 



similar to those of anul0-l, including leaves with promi- 
nent marginal teeth and paler rosettes, when compared 
with their wild-type Col-0 (Supplementary Fig. S3A-D at 
JXB online). In general, these phenotypes were less severe 
than those of anul0-l and were more apparent in younger 
rosettes, becoming indistinguishable from Col-0 14-15 das. 
The mutant phenotype of the Fj progeny from the crosses 
of anul0-l to SAIL_708_F05 and SAIL_659_F07 indicated 
that both T-DNA lines carry anulO alleles (Supplementary 
Fig. S3E, F), which were named anulO-2 and anulO-3 
(Fig. 3B). Because the anulO-2 and anulO-3 mutants carry 
T-DNA insertions in exons 3 and 1, respectively, which are 
expected to disrupt the coding potential of the mRNA, it 
is speculated that the milder effects of these mutant alleles 
might reflect differences between the Ler and Col-0 genetic 
backgrounds. 

To confirm further the identity of the ANU10 gene, a con- 
struct was made to express the coding region of Atlg28530 
constitutively. When transformed into anul0-l plants, the 
35S pm :ANU10 transgene fully complemented the defects in 
leaf shape (Fig. 1C, F), growth (Fig. 1C, J; Supplementary 
Fig. SI at JXB online), mesophyll development (Figs II, 2A, 
D, G, J, K), pigment levels (Fig. 1L), and H 2 0 2 accumulation 
(Supplementary Fig. S2) in each of six independent trans- 
formants. Therefore, the correct identification of ANU10 
as the Atlg28530 gene is supported both by the transgenic 
complementation studies and by the independent isolation of 
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mutations. Boxes and lines indicate exons and introns, respectively. 
A white box represents the 3'-UTR. Triangles indicate T-DNA insertions. 



three mutant alleles carrying lesions in the coding region of 
Atlg28530. 



The ANU1 0 gene encodes a protein of unknown 
function 

ANU10 is a single-copy gene in the nuclear genome of 
Arabidopsis. Because the gene is predicted to encode a protein 
of unknown function with no conserved domains, BLAST 
searches were carried out to identify similar protein sequences 
in public databases. Significant hits were found in the genomes 
of other higher plants and the moss Physcomitrella patens, 
but not in those of animals or other eukaryotes, including 
algae, suggesting that ANU10 belongs to a family of embry- 
ophyte-specific proteins (Supplementary Fig. S4 available at 
JXB online). The chloroplast localization of ANU10 was 
consistently predicted by several computational tools (see the 
Materials and methods), including TargetP (score=0.929) and 
Multiloc2 (score=0.57). ChloroP 1.1 predicted a chloroplast 
transit peptide in ANU10 and in most of its orthologues from 
other land plants (Supplementary Fig. S4A, Supplementary 
Table S2). In addition to the transit peptide, ANU10 and 
most of its orthologues were also predicted to have a trans- 
membrane domain (Supplementary Fig. S4A, Supplementary 
Table S2), suggesting that these proteins are anchored to chlo- 
roplast membranes. As an example, SOSUI predicted a trans- 
membrane domain spanning residues 421^143 in ANU10. 

To gain insight into the evolutionary origin of this protein 
family, HMMER searches were carried out using a profile 
made with the sequences of several ANU10 homologues from 
land plants. HMMER allowed the identification of some dis- 
tantly related sequences in Cyanobacteria. In line with these 
results, a search for known domains in the Pfam database 
(Punta et al, 2012) yielded a low significance hit to a domain 
of unknown function (DUF4335) present in some cyanobac- 
terial proteins. Together, these data indicate that ANU10 is 
conserved among land plants. However, unlike proteins such 
as CURT1A (Armbruster et al, 2013), which is functionally 
conserved in Cyanobacteria, the search for cyanobacterial 
orthologues of ANU10 did not yield obvious candidates. 



Analysis of ANU 10 expression 

To determine the expression pattern of ANU10 in wild- 
type plants, an ANU10 pro :GUS reporter transgene was 
constructed. Four independent transgenic plants stably 
expressing the transgene and showing the same staining 
pattern were studied. The GUS signal was broad in seed- 
lings collected 3 das (Fig. 4A), being particularly intense 
in incipient leaves (Fig. 4B). In plants collected 13 das, 
GUS activity was detected in roots, cotyledons, and leaves 
(Fig. 4C-H). In older plants, GUS activity was observed 
in stems, cauline leaves, flowers, and siliques (Fig. 4I-N). 
Remarkably, the GUS signal was more intense in young 
cauline leaves (Fig. 4J), flowers (Fig. 4K), and siliques 
(Fig. 4M) than in the corresponding mature organs 
(Fig. 41, L, N). The spatio-temporal expression analysis 
suggests that ANU10 is expressed in all plant organs and 
that its expression is particularly important in developing 
organs. 



ANU 10 localizes to plastids 

To determine experimentally the subcellular localization 
of the ANU 10 protein, an in-frame translational fusion of 
GFP to the C-terminal end of ANU 10 was made. Transgenic 
anul0-l plants expressing 35S pm :ANU10:GFP were used to 
visualize the GFP signal by confocal laser scanning micros- 
copy. The GFP signal was specifically detected in chloroplasts 
from four independent lines carrying the 35S pm :ANU10:GFP 
transgene (Fig. 5A-F). To exclude that the observed localiza- 
tion pattern represents an artefact due to overexpression of 
the ANU10:GFP fusion protein, the same translational fusion 
was also placed under the control of the endogenous ANU10 
promoter {ANU10 pro :ANU10:GFP). An identical distribution 
of GFP signal was observed in anul0-l plants expressing the 
ANU10 pro :ANU10:GFP transgene (Supplementary Fig. S5 at 
JXB online), demonstrating that the promoter chosen does 
not affect the subcellular localization of the fusion protein. 
Interestingly, both transgenes were able partially to rescue the 
anul0-l phenotype in young rosettes (Supplementary Fig. 
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Fig. 4. Visualization of ANU10 pro :GUS activity in a wild-type background. (A) Seedling, (B) detail of the shoot apex, (C) roots, and (D) cotyledon. (E-H) 
Vegetative leaves from the (E) first, (F) third, (G) fifth, and (H) seventh nodes. (I) Inflorescence, (J) young cauline leaf, (K) immature flowers, (L) mature 
flower, and (M) immature and (N) mature siliques. Pictures were taken (A) 3, (B) 6, (C-H) 13, and (l-N) 42 das. Scale bars indicate (A, D-G, I, J, L-N) 
1 mm, (C, H, K) 500 urn, and (B) 100 \xm. 



S6). The rescue was complete 21 das regardless of the pro- 
moter used, indicating that the ANU10:GFP fusion protein 
retains sufficient activity to complement the mutant pheno- 
type at this stage. To determine the suborganellar localization 
of the ANU10 protein, stromal and thylakoidal fractions were 
isolated from the chloroplasts of anulO-1 35S pm :ANU10: GFP 
transgenic plants. The presence of the ANU10:GFP fusion 
protein in each fraction was tested by western blotting using 
an anti-GFP antibody. The fusion protein was specifically 
detected in thylakoids, but not in the soluble stromal fraction 
(Fig. 5G), indicating that ANU10 is associated with thylakoi- 
dal membranes. 

Because ANU10 was found to be expressed in roots 
(Fig. 4C), root tissues were also examined under a confocal 
microscope. A strong GFP signal was detected in root tips from 
four independent lines expressing the ANU10 pro :ANU10: GFP 
transgene (Fig. 6A, B). The pattern of GFP signal paralleled 
the distribution of amyloplasts in root tips stained with lugol 
(Fig. 6C). Although amyloplasts play a central role in root 
graviperception, no obvious defects were detected in the grav- 
itropic response of anulO roots (Supplementary Fig. S7A 
at JXB online). Starch content was similar between anulO 
mutants and the wild type, as visualized after lugol stain- 
ing (Supplementary Fig. S7B). The primary root length 
was reduced in anulO-1 (1.75 ±0. 25cm) compared with her 
(2.15 ±0.25 cm) (Fig. 6D). Roots of anulO-2 and anulO-3 
were also shorter (2. 14 ±0.26 and 1.92 ±0.37 cm, respectively) 



than those of Col-0 (2.79 ±0.30 cm). However, this pheno- 
type, unlike the leaf phenotypes, was not restored in anulO-1 
35S pm :ANU10 transgenic plants, in which the mean length of 
the main root was 1.93 ±0.49 cm (Fig. 6D). 

Thylakoid membranes and LHCII trimers are less 
abundant in anul 0-1 

To investigate whether ANU10 deficiency affects the subor- 
ganellar organization of chloroplasts, the ultrastructure of 
chloroplasts from her, anulO-1, and tmulO-1 35S pro :ANU10 
mesophyll cells was studied using transmission electron 
microscopy. Chloroplasts were smaller in anulO-1 plants 
than in the wild type. While Ler chloroplasts are typi- 
cally lens shaped (Fig. 7A, B), anulO-1 chloroplasts were 
smaller and abnormally shaped (Fig. 7D, E). Fewer thy- 
lakoidal membranes were observed in anulO-1 chloroplasts 
(Fig. 7B, E). When observed at high magnification, it was 
found that anulO-1 thylakoids failed to stack and form 
typical grana. Under the conditions used here, wild-type 
grana were composed of 8.37 ±2.36 layers of thylakoidal 
membranes, while the number of layers per grana was only 
3.08 ± 1.02 for anulO-1 (Fig. 7F, L), as also reflected by the 
decreased height of anulO-1 grana (Fig. 7K). Moreover, the 
diameter of the grana and the number of grana per chloro- 
plast section were found to be reduced in anulO-1 (Fig. 7J, 
M). Chloroplast size, morphology, thylakoid abundance, 
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Fig. 5. Subcellular and suborganellar localization of ANU1 0. (A-F) 
Confocal micrographs of the subepidermal layer of palisade mesophyll 
cells from (A-C) Ler and (D-F) anul0-l 35S pro :ANU10:GFP transgenic 
plants. Micrographs show (A, D) the chlorophyll autofluorescence of the 
chloroplasts, (B, E) the GFP fluorescence, and (C, F) an overlay of the 
chlorophyll and GFP signals, showing their co-localization in (F). Pictures 
were taken from first-node leaves collected 1 6 das. Scale bars indicate 
50 |im. (G) Western blot analysis of the proteins in chloroplast, stroma, and 
thylakoid fractions isolated from anu10-1 35S pm :ANU10:GFP transgenic 
plants collected 1 6 das. Primary antibodies against GFP, the large Rubisco 
subunit (RbcL), and the PsbO subunit of PSII were used. Molecular mass 
markers are indicated on the right. 
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Fig. 6. Localization of the ANU1 0 protein and effect of anulO mutations 
on root growth. (A, B) Confocal micrographs of the root apex from 
an anu10-1 ANUW pm :ANU10:GFP transgenic plant: (A) differential 
interference contrast (DIC) image and (B) an overlay of the GFP 
fluorescence and the DIC image. (C) Root apex from a Ler plant stained 
with lugol. Scale bars indicate 30 urn. (D) Main root length ofanu10-1, 
anu10-2, and anu10-3 mutants, their respective wild types, and anu10-1 
35S pm :ANU10 transgenic plants. Error bars indicate standard deviations. 
Asterisks indicate values significantly different from the wild type in a f-test 
(P<0.05, n=30). 



and grana morphology and stacking were totally restored 
to those of the wild type in anulO-1 plants expressing the 
35S pm :ANU10 transgene (Fig. 7G-I, J-M). Because LHCII 
trimers are thought to participate in thylakoid stacking, 
blue native PAGE was used to study the levels of trimeric 
LHCII complexes in isolated thylakoids from 16 das plants. 
LHCII trimers migrated slightly more slowly and their levels 
were reduced in thylakoids of anulO-1 compared with those 
of her (Fig. 8A, lanes 1 and 2). The anulO-2 and anulO- 
3 mutants showed a similar reduction in the amount and 
mobility of LHCII trimers relative to the Col-0 wild type 
(Fig. 8A, lanes 3-5). The mobility and levels of LHCII 
trimers were fully restored in thylakoids of anulO-1 plants 
carrying the 35S pro :ANU10 transgene (Fig. 8A, lane 6). The 
levels of LHCII trimers were partially restored to wild-type 
levels in anulO-1 35S pm :ANU10:GFP plants (Fig. 8 A, lane 
7), in line with the partial phenotypic rescue seen in anulO-1 
35S pro :ANU10:GFP plants collected 16 das (Supplementary 
Fig. S6B at JXB online). Taken together, these results sug- 
gest that ANU10 is important for thylakoid biogenesis 
and grana formation, and do not allow the exclusion that 
a defect in the levels or composition of the LHCII trimers 
might be responsible for the reduced thylakoid stacking seen 
in the anulO-1 mutant. 



Expression of nuclear and chloroplast genes in 
anu10-1 

Quantitative RT-PCR (qRT-PCR) was used to explore 
whether ANU10 deficiency affects the expression of nuclear 
and plastid genes. The relative expression of nuclear genes 
encoding subunits of LHCII trimers (LHCB1, LHCB2 and 
LHCB3) was found to be ~0. 5-fold decreased in anulO-1 com- 
pared with her (Fig. 9A). Because the LHCB5 gene is known 
to be overexpressed in plants that do not express the LHCB1 
and LHCB2 proteins (Ruban et al, 2003), the expression of 
LHCB5 was also analysed. However, the expression levels of 
LHCB5 were not significantly different in anulO-1 and her 
(Fig. 9A). Because the impact of the anulO-1 mutation on 
the expression of LHCB genes might be an indirect conse- 
quence of a more general problem that triggers retrograde 
signals from the chloroplast to the nucleus, the expression of 
the HEMA1 gene, which is involved in chlorophyll biosyn- 
thesis and is particularly tightly regulated by chloroplast-to- 
nucleus signals (Mochizuki et al, 2008), was also examined. 
The expression of HEMA1 was significantly although mildly 
reduced in anulO-1 (Fig. 9A), suggesting that the reduced 
expression of LHCB genes is at least in part an indirect effect 
of the mutation on these signalling pathways. 



ANU10 is required for leaf development | 2399 




Fig. 7. Ultrastructure of anulO-1 chloroplasts. Transmission electron micrographs of palisade mesophyll cell chloroplasts from (A-C) Ler, (D-F) anu10-1 , 
and (G-l) anuW-1 35S pro :ANU10. Arrowheads in (F) indicate unstacked thylakoid membranes in the anuW-1 mutant. Pictures were taken from first-node 
leaves collected 16 das. Scale bars indicate (A, D, G) 5 urn, (B, H) 2 |im, (E) 1 (im, and (C, F, I) 200 nm. (J-M) Comparison of (J) the diameter (x-axis) and 
(K) height (y-axis) of granal stacks, (L) the number of membrane layers in granal stacks, and (M) the number of stacks in chloroplast sections from Ler, 
anu10-1 , and anu10-1 35S pm :ANU10. The number of grana per chloroplast section was determined from images of individual chroroplasts similar to 
those shown in B, E, and H. Asterisks indicate values significantly different from the wild type in a f-test (P<0.05, n=10-20). 



Because the phenotype of anulO-1 mutants might also be 
interpreted as a consequence of premature senescence, which 
might lead to grana disassembly and thylakoid swelling (Evans 
et al, 2010; Krupinska et al, 2012), the expression level of 



the nucleus-encoded ORESARA1 ( ORE1) transcription fac- 
tor, a senescence marker that has recently been reported to 
be a key regulator of the senescence response (Kim et al, 
2009; Balazadeh et al, 2010), was also monitored. ORE1 was 
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Fig. 8. Thylakoidal protein complexes in the anu10 mutants. Blue native PAGE of photosynthetic protein complexes from anu10-1 , anu10-2, and 
anuW-3 mutants, their respective wild types, and anuW-1 35S pro :ANU10 and anu10-1 35S pm :ANU10:GFP transgenic plants. PSII, photosystem II 
photosystem I; Cyt b B f, cytochrome b e f complex; LHCII, light-harvesting chlorophyll a/fa-binding protein complex II. 




expressed at similar levels in anulO-1 mutants and wild-type 
rosettes (Fig. 9A), suggesting that premature senescence is 
not the cause of the observed phenotype. 

To survey the effects of the anulO-1 mutation on the 
expression of plastid genes, the transcript levels of the plas- 
tid-encoded accD, rbcL, psbA, rrnl6, rrn23, and atpB genes 
were also quantified. A reduction only of rrnl6 expression in 
anulO-1 mutants was observed, while the expression levels of 
the other genes were not significantly different from those of 
Ler (Fig. 9B), indicating that the defects of anulO-1 chloro- 
plasts do not result from a general deregulation of organelle 
transcription. 



Discussion 

ANU1 0 encodes a novel plastid protein conserved 
throughout land plants 

The chloroplast proteome comprises ~1300 different proteins 
unambiguously identified in proteomic studies (Zybailov 
et al, 2008; Ferro et al, 2010), -92% of which are encoded 
in the nuclear genome. However, the function of a signifi- 
cant fraction (~30%) of these proteins remains unknown 
(Kleffmann et al, 2004), indicating that additional research 
is needed to understand fully the processes that take place 
in the chloroplast. In the present work, the ANU10 gene has 
been identified using a map-based cloning strategy. ANU10 
is the founding member of a plant-specific family of proteins 
that contain a chloroplast transit peptide and a transmem- 
brane domain. Using translational fusions to GFP, it has 
been shown that the ANU10 protein localizes to chloro- 
plasts. A transmembrane domain is also predicted to occur 
in most ANU10 homologues from other species. In line 
with this prediction, the ANU10 protein was detected using 



immunoblotting in the fraction corresponding to thylakoidal 
membranes, but not in the soluble (stromal) fraction derived 
from the chloroplasts. Study of the expression pattern and 
subcellular localization of the ANU10 protein shows that the 
gene is expressed in a variety of tissues, including cotyledons, 
leaves, flowers, and roots, and is consistent with microarray 
data available from the Arabidopsis eFP Browser (http:// 
bar.utoronto.ca/efp/cgi-bin/efpWeb.cgi). In close agreement 
with this expression pattern, GFP signal corresponding to 
the ANU10:GFP fusion protein was detected not only in 
the leaf chloroplasts, but also in the root amyloplasts. The 
data point to a differential requirement for ANU10 protein 
in roots and leaves. Although the expression of the full-length 
coding sequence of ANU10 restored the wild-type phenotype 
in anulO-1 mutant leaves, the root length was not restored 
to wild-type values, suggesting that roots are sensitive to an 
excess of ANU10 protein. 

ANU1 0 is required for thylakoid biogenesis and grana 
stacking 

Very little is known as yet about the mechanisms that lead 
to normal thylakoid biogenesis and grana stacking. The pre- 
sent finding that anulO-1 mutants exhibit aberrant thylakoid 
stacking and lack the typical grana seen in normal chloro- 
plasts is expected to help in furthering understanding of these 
processes. Although a mechanistic link between the molecular 
function and the observed phenotypes is still lacking, ANU10 
might contribute to the characteristic internal organization 
of chloroplasts in several different ways. 

First, previous authors have proposed that LHCII trimers 
play a crucial role in thylakoid stacking and grana formation 
(Day et al, 1984; Allen and Forsberg, 2001; Cui et al, 2011). 
A lower amount of LHCII trimers in anulO-1 thylakoids was 
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Fig. 9. Expression of nuclear and plastid genes in the anu10-1 mutant. 
(A, B) qRT-PCR analysis of the expression of (A) LHCB1 , LHCB2, LHCB3, 
LHCB5, HEMA1, and ORE1 nuclear genes, and (B) accD, psbA, rbcL, 
rrnW, rrn23, and atpB plastid genes in Ler and anuW-1 rosettes collected 
1 6 das. Bars indicate the relative expression levels, determined by the 
comparative C T method, and normalized with the expression of the 18S 
rRNA housekeeping gene. Error bars indicate the interval delimited by 
2" (AACT±SD| . Asterisks indicate ACj values significantly different from those 
of Ler in a Mann-Whitney U-test (P<0.01 ; n=9). 

detected here, in line with the proposed role of trimeric LHCII 
complexes in thylakoid stacking. Reduced levels of LHCII 
trimers in the gdcl-3 mutant have been linked to impaired 
grana stacking in Arabidopsis (Cui et al, 2011), and other 
authors have shown that normal thylakoid stacking can still 
occur in plants with reduced levels of LHCB1 and LHCB2 
(Andersson et al, 2003), which are compensated by higher 
levels of LHCB5 (Ruban et al, 2003). Stacks of thylakoidal 
membranes have been found in other mutants with reduced 
levels of LHCB transcripts, such as genomes uncoupled! 
(gunl) (Susek et al , 1993), immutans (Wetzel et al. , 1994), and 
CAB underexpressed4 (cue4) and cue9 ( Lopez -Juez et al, 
1998). Therefore, it is presently unclear if the moderate reduc- 
tion in LHCB transcription and LHCII trimers detected in 
anulO mutants can account for the thylakoid phenotype. 

Secondly, the lack of grana stacks in anulO mutants might 
be a consequence of a more general deficiency in thylakoid 
biogenesis. Thylakoidal membranes, which are less abundant 
in anul0-l chloroplasts, are thought to originate from vesicles 



that form at the inner envelope of chloroplasts and subse- 
quently incorporate into developing thylakoids (Eggink et al. , 
2001; Tanz et al, 2012). Three proteins known to participate 
in vesicle formation and transport during thylakoid biogen- 
esis are VESICLE-INDUCING PROTEIN IN PLASTIDS1 
(VIPP1), CHLOROPLAST SECRETION-ASSOCIATED 
RAS1 (CPSAR1), and SNOWY COTYLEDON2/SHI- 
YO-U1 (SC02/CY01). The VIPP1 protein localizes to 
the thylakoidal and chloroplast inner membranes, and the 
high-chlorophyll fluorescence!55 (hcfl55) mutant shows that 
VIPP1 is required for the formation of the vesicles (Kroll 
et al, 2001). CPSAR1 dually localizes to the chloroplast 
stroma and the inner envelope, and has been proposed to 
participate in the initiation of vesicles from the inner enve- 
lope (Garcia et al, 2010). SC02/CY01 is required for the 
trafficking of vesicles from the inner envelope to developing 
thylakoids (Tanz et al., 2012). Interestingly, the anul0-l phe- 
notype resembles the phenotype of loss-of-function hcfl55, 
cpsarl, and sco2-l mutants as regards the low abundance of 
thylakoids (Kroll et al, 2001; Garcia et al., 2010; Tanz et al, 
2012). This raises the possibility that ANU10 participates in 
thylakoid biogenesis by modulating vesicle integration into 
developing thylakoids. Under this scenario, the reduced lev- 
els of LHCII trimeric forms detected in anul0-l thylakoids 
would simply be an indirect consequence of inefficient thy- 
lakoid biogenesis, as LHCB proteins have been proposed to 
be targeted to thylakoids, at least in part, through their incor- 
poration into developing vesicles in the inner membranes of 
chloroplasts (Eggink et al, 2001; Tanz et al, 2012). 

Because some proteins involved in plastid gene expression, 
such as pTAC14 (Gao et al, 2011) or the nucleus-encoded 
RPOTmp polymerase (Azevedo et al, 2008), are associated 
with the thylakoidal membranes, and ptacl4 mutants have 
been reported to lack grana and thylakoidal membranes 
(Gao et al. , 201 1), it was also considered here that the defects 
in thylakoid abundance and stacking might result from a 
more general problem in plastid gene expression in anulO- 
1. However, the results of the present study indicate that the 
expression of plastid genes is not generally affected. 

Leaf senescence has also been associated with thyla- 
koid defects in higher plants (Evans et al, 2010; Krupinska 
et al, 2012). In the anul0-l mutant, however, the senescence 
marker ORE1 was expressed at normal levels, suggesting that 
senescence is not the primary cause of the observed plastid 
phenotypes. 

Loss ofANUW function compromises leaf 
development 

Altered leaf anatomy is a trait common to many mutants car- 
rying lesions in nuclear genes that encode plastid proteins. 
Some examples are the immutans (Wetzel et al, 1994; Carol 
et al. , 1999; Aluru et al. , 2001), cuel (Li et al , 1995; Streatfield 
et al, 1999), pale cressl {pad) (Reiter et al, 1994), and 
rugosa2 (rug!) (Quesada et al. , 201 1) mutants of Arabidopsis, 
the differentiation and greening (dag) mutant of Antirrhinum 
majus (Chatterjee et al, 1996), and the defective chloroplasts 
and leaves (del) mutant of tomato (Keddie et al , 1996). These 
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mutants display variegated or pale-green phenotypes and 
have altered plastid and mesophyll development. In some var- 
iegated mutants, such as rug2, the pale sectors exhibit severe 
developmental defects in plastids and mesophyll cells, while 
the cells and chloroplasts within green sectors are usually less 
affected (Quesada et al, 2011). 

Because some mutations damaging plastid-localized pro- 
teins affect both plastid development and mesophyll cell dif- 
ferentiation, plastid development and leaf morphogenesis 
have been hypothesized to be tightly coordinated processes 
via retrograde plastid-to-nucleus signalling (Rodermel, 2001). 
The signal is assumed to originate from chloroplasts with 
arrested differentiation or altered photosynthetic metabo- 
lism, and modulates the expression of nuclear genes encoding 
plastidial and other proteins (Susek et al, 1993; Koussevitzky 
et al, 2007). Pyke et al (2000) proposed that a cell-autono- 
mous signal originating in the plastids is required for normal 
development of palisade mesophyll cells. Like the above- 
mentioned mutants, the anul0-l mutant exhibits altered mes- 
ophyll cell development, with larger and less densely packed 
palisade cells. The present morphometric analysis shows that 
anul0-l leaves have fewer, larger mesophyll cells per area unit, 
suggesting either a premature transition from proliferation to 
cell expansion, or a defective cell proliferation. The increased 
extent of expansion observed in anul0-l is in contrast to pre- 
vious studies (Andriankaja et al, 2012), which have shown 
a delay in mesophyll cell differentiation when chloroplasts 
are prevented from differentiating. Different chloroplast-to- 
nucleus signals acting in differentiating and mature chloro- 
plasts might explain the different effect on mesophyll cell 
expansion. 

LHCB genes, among many other nuclear genes related to 
plastid photosynthetic machinery and other diverse functions, 
are known to be differentially expressed in response to chlo- 
roplast retrograde signalling (Susek et al, 1993; Koussevitzky 
et al, 2007). Furthermore, recent experimental evidence sug- 
gests that retrograde signalling also modulates the expression 
and activation of cyclin-dependent kinases, which are key reg- 
ulators of the cell cycle (Kobayashi et al, 2009; Andriankaja 
et al. , 2012). The defects of anul0-l chloroplasts might trigger 
a retrograde signal leading to changes in nuclear gene expres- 
sion, as exemplified by the reduced expression of HEMA1 
and the LHCB1, LHCB2, and LHCB3 genes. 

Concluding remarks 

A nuclear gene, ANU10, has been identified, whose loss-of- 
function mutations lead to pleiotropic defects in leaf develop- 
ment and plastid internal organization, including a dramatic 
reduction in thylakoidal membranes and grana stacking. The 
results indicate that these phenotypes are not merely a con- 
sequence of a premature senescence response or of a general 
defect in the transcription of plastid-encoded genes. The 
expression of several nuclear genes, including HEMA1 and 
several LHCB genes, was found to be altered in the anulO- 
1 mutant, suggesting that the plastid defects trigger a retro- 
grade (chloroplast-to-nucleus) signal that might account for 
the mesophyll phenotype. 



Supplementary data 

Supplementary data are available at JXB online. 
Figure SI. Growth of the anul0-l mutant. 
Figure S2. ROS in amtl0-l leaves. 

Figure S3. Rosette phenotype of homozygotes and het- 
erozygotes for the anulO alleles. 

Figure S4. Phylogenetic analysis of ANU10 and its puta- 
tive orthologues. 

Figure S5. Subcellular localization of the ANU10:GFP 
fusion protein in the anul0-l ANU10 pw :ANU10:GFP 
transgenic line. 

Figure S6. Phenotypic rescue of anul0-l by the 
ANU10:GFP fusion protein. 

Figure S7. Growth and starch content of anulO roots. 

Table SI. Primer sets used in this work. 

Table S2. Transmembrane domains and chloroplast 
transit peptides in the ANU10 protein and its putative 
orthologues. 
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